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We investigate Aharonov-Bohm osciUations of the current 
through a strongly correlated quantum dot embedded in an 
arbitrary scattering geometry. Resonant-tunneling processes 
lead to a flux-dependent renormalization of the dot level. As 
a consequence we obtain a fine structure of the current oscil- 
lations which is controlled by quantum fluctuations. Strong 
Coulomb repulsion leads to a continuous bias voltage depen- 
dent phase shift and, in the nonlinear response regime, de- 
stroys the symmetry of the differential conductance under a 
sign change of the external flux. 



Phase-sensitive transport properties of interacting 
mesoscopic systems are important for several reasons. 
The small size of the samples gives rise to capacitances 
of order 1Q~^^F which induce Coulomb blockade effects 
101 and demand the necessity to generalize the Landauer- 
Biittiker formalism to systems with strong interac- 
tions. Furthermore, the investigation of Aharonov-Bohm 
oscillations through quantum dots with strong Coulomb 
repulsion might give further experimental evidence for 
resonant tunneling and Kondo phenomena in nonequi- 
librium systems 

Interference effects in the Coulomb blockade regime 
have been measured by Yacoby et al. ||] by studying a 
quantum dot embedded in an Aharonov-Bohm ring. This 
experiment demonstrates that phase-coherent transport 
through quantum dots is possible in realistic experiments 
and is not destroyed by inelastic interactions. Recent the- 
oretical work on Aharonov-Bohm oscillations in a meso- 
scopic ring with a quantum dot uses a noninteract- 
ing model. Using the symmetry of the current under 
sign change of the external flux in the linear response 
regime HJl^ it was shown in Ref. that the phase of 
the Aharonov-Bohm oscillations can only take two pos- 
sible values as function of the gate voltage on the dot. 
However, the experiment of Yacoby et al. is performed 
in the Coulomb blockade regime where interaction effects 
are important. In this letter we will take such correla- 
tions into account by setting up a complete and general 
theory for interference phenomena in strongly interacting 
quantum dots embedded in an arbitrary noninteracting 
multi-probe and multi-channel scattering structure. As a 
consequence we will show that the symmetry under sign 



change of the external flux is broken in the nonlinear 
response regime and that the phase can change continu- 
ously as a function of the bias voltage. Furthermore, we 
will analyze in detail the current oscillations as a func- 
tion of the gate voltage caused by the flux-modulated 
renormalization of the local energy level of the dot. 

To have a specific example, we will study the system 
shown in Fig. ^ although our formalism is valid for an ar- 
bitrary scattering geometry. For simplicity we start with 
the case of one-channel leads. The system without the 
quantum dot is described by scattering waves with zero 
boundary conditions at the tunneling barriers of the dot. 
Thus, in energy representation, this part of the Hamilto- 
nian is given by Hs = J2aa I * <^ "Ll'^) aQ^le) , where 
ajjg.(e) creates an incoming scattering wave in probe a 
with spin cr and total energy e. The isolated dot is de- 
scribed by Hd = J2a ^(y^ld^ + UY,a«j' "cr"cr' ^ith 

single particle energies and on-site repulsion U . The 
position of the dot levels are controlled by an external 
gate voltage and U ^ 1 — 5K corresponds to the charging 
energy 0. 

The tunneling of the electrons into or out of the dot is 
described by 



Ht^^ de {ia(e)aL(e)d, + /i.e.} 



(1) 



Here, ta{e) — X]i=L/i? ("^^l^^) ^'^s ^^le tunneling ma- 
trix elements in energy representation, where ti are real 
quantities and {x\ae) is the spin-independent scatter- 
ing wave from reservoir a with energy e at position x. 
By Xi , i = L, R, we denote an arbitrary point in the 
one-dimensional left or right lead which is connected 
to the dot Due to zero boundary conditions we 

have {xi\ae) — p(e)^/^^^(e) sin(fc(e)a:i) with the one- 
dimensional density of states p(e) — l/{TThv{e)) and en- 
ergy e = h^k{ef/{2m) = \mv{ef . The coefHcients 
depend on the detailed scattering problem under consid- 
eration. We have chosen the tunneling matrix elements 
ti as real parameters which means that we shift the com- 
plete flux dependence to the scattering Hamiltonian Hs 
via a standard gauge transformation. 

Following Biittiker we will use the following rep- 
resentation of the current operator in probe a 
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where 6acr(e) = ■Sa/3(e)ci/3CT(e) annihilates an outgoing 
carrier in probe a and s is the scattering matrix of the 
system without the dot. To calculate the average current 
la = {la) in the stationary limit we need the station- 
ary real-time Green's function G<{E) = J dte'^*G<(t) 
in Fourier space of two scattering field operators: 
Ga<T,a'a'(£'e';^) = «(aL(e>^)aa'<T'(e')>- Using the matrix 

notation G = ( ^ j , where G^ and G^ are the 

retarded and advanced Green's functions, and applying 
the Keldysh technique we can express the Green's 
function Gaa,a'a cxactly by the local Green's function 
Ga of the dot, Gaa,a'a{<^, \ E) = ga{e; E)6a,a'S{e - e') + 
ta{e)ta'{e-')*ga{^\E)Ga{E)ga'{e.'\E) , where we have al- 
ready used spin conservation. The Green's functions 
ija correspond to the Hamiltonian Hs and are given 
by g^'^{e-E) = [E - e ± zO+)-i and g<{e-E) = 
'^'^ifa{E)d{E^e) where fa is the Fermi distribution func- 
tion of reservoir a. Using this result in calculating the 
average current, inserting the form of the tunneling ma- 
trix elements and performing the energy integrations , 
we obtain 



(3) 



where la'^ is the current without the dot (given by 
the usual Landauer-Biittiker formula) and Aaa' = 
E„-(r.U,)i^r^i' with r,(e) = 2^p(e)tf . 

Eq. (ph is the first central result of this paper. It re- 
lates the current in probe a exactly to the local Green's 
fimctions of the dot and the scattering properties of the 
noninteracting medium surrounding the dot. This for- 
mula is a natural generalization of the Landauer-Biittiker 
formula to an interacting quantum dot. Furthermore it 
generalizes current formulas through quantum dots con- 
nected to two leads without any possibility of a direct 
transition between the probes Here we are able to 
account for such transitions opening the possibility to 
study interference phenomena in the presence of locally 
interacting subsystems. The generalization of Eq. (^) to 
multi-channel leads is straightforward. Again following 
Ref. [|l3| the field operators Uaa- in Eq. (||) have to be 
treated as vectors with a channel index n. Equivalently, 
the matrix elements Saf3 and Aa/3 have to be treated like 
Za X Zfj matrices where Za is the number of transverse 
channels in lead a. The final formula for the current is 
then exactly like Eq. (||) except that we have to take the 
trace of the matrix multiplication slaSa-yA^p. 

The scattering matrices in Eq. (0) can be found by 
straightforward quantum-mechanical considerations de- 
pending on the specific geometry. For the Green's func- 
tions of the dot we will use a real-time technique devel- 
oped in Ref. which has been applied to a quantum 



dot in Ref. For a degenerate dot level (i.e., = 
independent of spin) and in the U = oo limit, one obtains 

{E) ^2T:i"f^{E)\E - ea- cr{E)\-^. Here, 
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a{E) = / dE 



^,M ^+{E')+^-{E') 
E-E' + iQ+ 



(4) 



has the form of a self-energy which describes the renor- 
malization and broadening of the dot level e^, 7*(£') = 
Y.a\UE)?ft{E) = l/i^7r)j:^AaaiE)f^iE) is the 
classical rate for a particle tunneling in or out of the dot, 
and /+ = fa while /~ = 1 — fa- The retarded Green's 
function follows from Im G" ^ l/(2i ) (G> - G< ) and the 
real part is obtained from the Kramers-Kronig relation. 

The explicit result for the Green's functions together 
with the expression (^ for the current constitutes a com- 
plete theory of interference effects in mesoscopic scatter- 
ing geometries with an interacting part given by a quan- 
tum dot with one degenerate level. Our result satisfies 
current conservation J^a ~ 0, and all currents vanish 
in equilibrium. Furthermore, for the special case M = 1 
where the Coulomb interaction does not play any role, 
our result is exact and can be shown to agree with the 
Landauer-Biittiker formalism. 

The real part of the self-energy a describes the renor- 
malization of the dot level. If we neglect the energy 
dependence of Aaa at the Fermi level, we obtain from 
Eq. (^) for a two-terminal system 

M -1 , , 

Rea^Recri + — [{An + A22 Xi + X2 + 

Stt 

+ {All - A22){xi - X2)] , (5) 

where <ti is the self-energy for M = I and XaiE) = 
Re J dE'fa{E')/{E-E'+iO+). Using a Lorentzian cutoff 
at D (which will be of the order of the Coulomb repulsion 
[/), we obtain Xa{E) = V^(i + - Re 7/;(i + i^^) 
where ip is the digamma function and fia the chemical 
potential of reservoir a. oi is always a symmetric func- 
tion of the external flux <&. Furthermore, for a spatially 
symmetric situation as in Fig. ^ An ± is an even 
(odd) function of the phase ~ 27r$/$o (^0 being the 
flux quantum). Due to Re ai the level position of the 
dot will oscillate with tp with an amplitude of the order 
of r and phase or tt. For M > 1, there can be loga- 
rithmic corrections in temperature and bias voltage for 
the amplitude and phase of this oscillation due to the Xa 
functions. The latter terms usually lead to Kondo-like 
correlations [^[^]. 

To exhibit the consequences of the oscillation of the 
renormalized dot level we will now apply our results to 
the specific scattering geometry of Fig. ^ which corre- 
sponds to the experimental setup of Ref. For sim- 
plicity we assume a one-dimensional structure and we 
use the same scattering matrices s*'° for the incoming 
and outgoing junctions as in Ref. joj. The scattering 
matrix of the upper arm (including the flux and the 
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phases accumulated by free motion) is written in the form 
= p ^ ^^i^ ^ ^ , where p = e'*^' is the phase ac- 
quired by free motion through the upper arm. Further- 
more, we take the length of the leads connected to the 
quantum dot a,s II = Ir = and we assume a symmetric 
quantum dot with = Tji = T. 

We will look explicitly at two cases, viz., perfect trans- 
mission through the upper arm given by r = 0, t = 1, or 
weak transmission described hy r = ~l, t = i\t\. In the 
first case we obtain after a straightforward calculation 

Sll = S22 = ^Pip- 1), Si2{(f') = S2li-(p) = \p{p+ l)e"^, 

A{ = AJ^ = \i^{2 - p) and A{{^) = (-^) = 
— ■^ipy/pe^'^ . In the second case one obtains sn = S22 — 
-p, si2((y3) = S2i(-<p) = 5^p|^|e^'^, = A§ = 
and A^iip) = A^{~ip) = ^j^\t\e"f. The phase p can- 
not be determined and will have some specific value in 
the experiment. We assume here always to be in the 
quantum region, i.e., the lengths associated with tem- 
perature, bias voltage, and F should exceed the system 
length, so that we can neglect the energy dependence of 
p. In the perfect transmission case we choose p = i and 
get Re a\ = — |F(3-|-cos(/?), A11+A22 = T{3 + cosLp) and 
All — ^22 = — 2Fsin(/3. For weak transmission we take 
p — 1 and obtain Re cri — ~^\t\T cos ip, An + A22 = 2F 
and All — ^22 ~ — \t\T sin ip. 

In Fig. H we show the linear conductance for perfect 
transmission and T = F = 0.01 (in units of the cutoff D) 
for various positions of the dot level. We have assumed 
M — 2, i.e., the interacting case. The linear conductance 
is symmetric under a sign change of the flux since the 
sin (p term is absent in Eq. (||) for xi = X2- If the position 
of the dot level is below « —0.02, the current has a 
maximum around if = 0. For higher values of e^, the 
current has a minimum at (/? = 0. Although this looks 
similar to the abrupt phase change of it described in Ref. 

, Fig. 1^ shows that the response of the system cannot 
be described by the concept of a "phase shift". What 
happens instead is that the current (as a function of ip) 
changes its functional form. The scale of the transition is 
independent of temperature and is given by the intrinsic 
parameter F. For higher temperatures we find a smaller 
amplitude of the current oscillations but the qualitative 
picture remains. For weak transmission, the results are 
similar but the scale is given by tT. 

We will now turn to the nonlinear conductance. Figure 
^ depicts the differential conductance for weak transmis- 
sion as a function of (p for different voltages and level 
positions (T = lOF = 0.1). For = 0, the differential 
conductance is symmetric around ip ~ 0, in the nonlin- 
ear response case, this symmetry is absent due to the 
last term in Eq. (|^). The behavior of the differential 
conductance at the origin changes from a minimum to 
a maximum as a function of the level position e^, this 
time rather abruptly (energy scale tT). The asymmetry 



of the conductance curves for finite voltages is a genuine 
interaction effect; it disappears for M = 1. Furthermore 
we observe a continuous phase shift of the current os- 
cillations as function of the bias voltage which again is 
absent for M = 1. It is determined by the last term in 
Eq. (^ as well as by simp terms occuring explicitly in 
the current formula (||) via the A^^ matrices. Note that 
the temperature is one order of magnitude larger than 
F in this figure, i.e., an interference experiment of this 
type might yield information about correlation effects at 
temperatures which are accessible in experiments p8[ |. 

Finally we want to comment on the influence of in- 
teractions on the relative phase of the Aharonov-Bohm 
oscillations at successive peaks in the linear conductance 
as function of the gate voltage. In a noninteracting model 
two adjacent peaks correspond to transport through two 
different energy levels of the dot which have different par- 
ity. Thus the relative sign of and tji would change 
from one level to the next and consequently one ex- 
pects a phase shift of tt. However, in the experiment 
of Yacoby et al. no phase shift was measured. In ad- 
dition to the discussion of Ref. [Q, a strong Coulomb 
repulsion on the quantum dot could be an explanation 
for this observation. If there are N states on the dot 
which lie close together in energy but with the same par- 
ity in longitudinal direction (e.g. spin degenerate states 
or states differing in the transverse channel number), 
there would be N adjacent Coulomb peaks with the same 
phase of the Aharonov-Bohm oscillations. The distance 
of these Coulomb peaks is given by the charging energy U 
whereas in the noninteracting case all these peaks would 
fall together into one single peak. Therefore we conclude 
that in the presence of interactions the parity of the en- 
ergy levels contributing to transport at adjacent Coulomb 
peaks can be the same which provides an explanation for 
Yacoby's experiment. 

In conclusion, we have presented a complete theory for 
interference phenomena in strongly correlated quantum 
dots embedded in a scattering geometry. On one hand, 
we have found that the functional form of ^ {(p) is chang- 
ing with the gate voltage on the scale of temperature- 
independent intrinsic parameters. In linear response this 
change cannot be interpreted as a phase shift. On the 
other hand, we have shown that in the nonlinear response 
regime, correlation effects break the symmetry under sign 
change of the external flux and lead to a real continuous 
phase shift as a function of the bias voltage. 
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FIG. 1. Geometry of the model system studied here. Leads 
1 and 2 connect to the left and right reservoir (shaded). The 
ring is connected to a quantum dot via high tunneling barri- 
ers. $ is the flux penetrating the ring. 
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FIG. 2. Linear conductance (in units of e^//i) as a function 
of magnetic flux for various positions of the level in the dot 
(T — r = 0.01, M — 2). We have assumed perfect transmis- 
sion through the ugper arm of the system as in the experiment 
by Yacoby et al. 



FIG. 3. Differential conductance (in units of e^ /h) for the 
case of weak transmission (f = 0.1) as a function of magnetic 
flux for various voltages and positions of the level in the dot 
(T = lOr = 0.1, M = 2). Note that the linear conductance 
(V^ = 0) is symmetric around = 0, whereas there is no such 
symmetry in the nonlinear response case. 
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